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Abstract 

We consider a one-dimensional loop of circumference L crossed by a constant magnetic flux $ 
and connected to an infinite lead with coupling parameter e. Assuming that the initial state 
ipQ of the particle is confined inside the loop and evolves freely, we analyse the time evolution 
of the nonescape probability P{ipQ,L,^,e,t), which is the probability that the particle will still 
be inside the loop at some later time t. In appropriate units, we found that P(tl'o,L,^,e,t) = 
Pooi'ipo,^) + XlfcLi CfcCV'O) -f') ^) The constant Poo{ipOi^) is independent of L and e, and 

vanishes if ipo has no bound state components or if | cos(<I>)| 7^ 1. The coefficients Cl{^l^o,L,^,£) 
and C^{ij)Q,L, e) depend on the initial state "00 of the particle, but only the momentum k = ^/L 
is involved. There are initial states ipo for which P{iIjq, L,^,e,t) ~ CsHjo, L,^,e)/t^ , as t — )• 00, 
where (5 = 1 if cos($) = 1 and 6 = 3 if cos(<I>) 7^ 1. Thus, by submitting the loop to an external 
magnetic flux, one may induce a radical change in the asymptotic decay rate of P{'tpQ, L,^,£,t). 
Interestingly, if cos(<I>) = 1, then Ci{ipo, L,^,e) decreases with e (i.e. the particle escapes faster in 
the long run) while in the case cos($) 7^ 1, the coefficient C■s{^lJo,L,^,e) increases with e (i.e. the 
particle escapes slower in the long run). Assuming the particle to be initially in a bound state of 
the loop with $ = 0, we compute explicit relations and present some numerical results showing a 
global picture in time of P{tpo, L, e, t). Finally, by using the pseudo-spectral method, we consider 
the interacting case with soft-core Coulomb potentials. 

PACS numbers: 03.65.-w, 03.65.Db, 03.65.Nk 
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I. INTRODUCTION 



The study of open quantum systems has attracted many researchers in theoretical physics 
during the last decades and is still under intense investigation. One of the most basic 
problems in this field is to understand how a quantum particle initially confined inside an 
open cavity will escape (see e.g. [ij] for a general discussion). In the literature, such systems 
are commonly described in terms of the following quantities: the nonescape probability P{t), 
which is the probability that the particle will still be inside the cavity at some later time t 
and the survival probability S(t), which is the probability to find the particle in the initial 
state at time t Let us suppose that the configuration space is the one-dimensional half- 
line [0,00) and that the cavity is the interval [0,L]. Then, at time t = 0, the particle is 
supposed to be in a normalized state ipo satisfying ipo{x) = for all x ^ [0,L]. If H denotes 
the Hamiltonian of the particle and ip{x,t) = {e~^^^^'^ijjQ){x) its state at time t > 0, where h 
is the reduced Planck constant, then the nonescape and survival probabilities are given by 



P{t) = I ijj{x,t)ijj{x,t)dx , (1) 
-L 2 







S{t) 



ipo{x)ilj{x, t)dx 







(2) 



where the bar " " denotes complex conjugation. It is easy to check that these quantities 
satisfy S{0) = P(0) = 1 and < S{t) < P{t) < 1 for all times t > 0. Ahhough both 
quantities are interesting, we shall discuss only P{t) in this paper. Note, however, that by 
analysing the time decay of P(t) one automatically gets an upper bound for S(t). 

The nonescape probability P(t) has been discussed in free ID quantum systems [sl, in 
free 2D quantum billiard systems ^ , [sj and also in 3D quantum systems with finite range 
potentials {g]. Interestingly, if the initial state ipo has no bound state components, then it is 
always found that the long-time behaviour of P{t) is a power law: P{t) = C/t^ + 0{l/t^~^^) 
as t — 00, where the constant C and the exponent S depend on the geometry of the system 
and on the initial state tpo of the particle. Some readers may be surprised since it is well 
known that various spontaneous decays have been observed to follow an exponential law. 
This apparent discrepancy may actually be only a time scale problem. Indeed, although 
quantum mechanics clearly predicts in general some deviation from exponential decay, the 
following phases may well occur ^-9\- for a short time the decay is parabolic P{t) = 1 — ct^, 
then it is approximately exponential and finally becomes a power-law at very large times. 



In the present paper, we consider the model introduced in [10| in which the cavity is 
a one-dimensional loop [0, L] crossed by a constant magnetic flux $ and connected, with 
coupling parameter e, to an inflnite one-dimensional lead [L, oo) (the points a; = and 
X = L are identifled), so that the particle may escape freely from the loop. See Fig. [1] This 
model contains rich physics and is simple enough to be amenable by analytical means or 
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. Nevertheless, 



numerical simulations, and thus has been the subject of several works 
the nonescape probability P{tljQ, L,^,e,t) has never been discussed in this model and it is 
interesting to know precisely how it depends on the initial state ipQ of the particle and on 
the three physical parameters: the loop's length L, the external magnetic flux $ and the 
coupling e. 




FIG. 1. A loop connected to a lead. 



In appropriate units, here is a summary of our main results. 

• Let A > 0, then one has the following scaling law: 

Piijo, AL, $, e, t) = Piijo. L, $, e, ^) . (3) 

Thus, increasing for example the length of the loop by a factor A = 2 leads to a 
nonescape probability evolving A^ = 4 times slower. 

• The nonescape probability may be written as 

P(^o, L, $, t) = PUi'o, $) + f: . (4) 

k=l 

The constant Poo('^0; ^) is given in (158|) : it is independent of L and e, and vanishes if 
ipo has no bound state components or if | cos(<I>)| ^ 1. In the inflnite series, the leading 
term always has an odd power. The coefficients Ci{ipQ, L, $, e) and C^lipQ, L, $, e) are 
given in ( l^ -f p5|l : they depend on the initial state iI^q of the particle, but only the 
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momentum /c = $/L is involved. When $ = 0, this coincides with the well-known 
fact that only the zero momentum of ipo plays a role in the long-time behaviour of 
the nonescape probability jsj. If cos($) = 1, then there are initial states tpo with 
no bound state components [e.g. (185|) ] for which Ci{ipo, L,^,e) 7^ 0, while the coeffi- 
cients Ci{tpQ, L,^,e) and C2{ipo, L,^,e) always vanish if cos($) 7^ 1. The important 
consequence is that there are initial states ipQ giving a power law decay: 

m..L.^.M)- ^'"^'°';-^'^' +0(^) . (5) 

where 5 = 1 if cos($) = 1 and 5 = 3 if cos(<l') 7^ 1. In particular, we see that by 
submitting the loop to an external magnetic flux (thus going from $ = to $ 7^ 0), 
one may induce a radical change in the decay rate of P{ipQ, L,^,e,t). If cos($) = 1, 
then Ci{ipo, L,^,e) decreases with e (i.e. the particle escapes faster in the long run) 
and if cos($) 7^ 1, then Cs{ipo, L, $, e) increases with e (i.e. the particle escapes slower 
in the long run). Interestingly, similar features were obtained in a circular dielectric 



cavity containing classical waves 



191], the nature of the waves (TM or TE) and the 



refractive index of the cavity playing a similar role than $ and e. 

• Assuming the particle to be initially in a bound state of the loop with $ = [see (172]) ]. 
we computed explicitly P^{ipo,^), Cii^po, L,^,e) and Csi^po, L,^,e) [see fl75]) - fl75]) ]. 
As one may expect, we found that higher energetic bound states decay faster. Interest- 
ingly, we observed that the coefficient C^lipo, L,^,e) oscillates, non-periodically, with 
$; see Fig. [51 We also present some numerical results showing a global picture in time 
of P{ipo, L,^,e,t); see Figs. [61-fIOl They reveal in particular that if cos($) 7^ 1, then 
by increasing the value of e, the particle escapes faster in the beginning but slower in 
the long run as stated previously. 

This paper is organized as follows. In Sec. [TTl we model the situation depicted in Fig. [T] in 



the framework of standard quantum mechanics. In the paper [10|, the scattering processes 
occurring at the connecting point (between the loop and the lead) are described in terms of 
an energy-independent scattering matrix S. We shall show that this scattering matrix may 
be used to prescribe boundary conditions at the connecting point for which the Hamiltonian 
of the particle is self-adjoint. In Sec. IIIIl we solve the time-dependent Schrodinger equation 
and derive its associated propagator. We then analyse the long time behaviour of the 
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propagator and consequently that of the nonescape probabihty. We discuss our general 
results in Sec. IIVI and consider some particular initial states in Sec. El In Sec. IVll we 
present an interesting and highly non-trivial application of the pseudo-spectral method by 
considering the interacting case with soft-core Coulomb potentials. Finally, in Sec. IVIII we 
make some concluding remarks. 



II. THE HAMILTONIAN 



In standard quantum mechanics, the situation depicted in Fig. [T]is modelled as follows: 
the physical states of the particle belong to the Hilbert space T-L = L'^{[0, L]) © L'^{[L, oo)) 
and the Hamiltonian of the particle reads 

H= p^-fx[o,L](Q)]' , (6) 

where X[o,l](Q) is the operator of multiplication by X[o,l]{x), where xio,l]{x) = 1 if x G [0, L] 
and X[o,L]ix) = otherwise. Here we have set h = 1 for the reduced Planck constant, e = 1 
for the electric charge, c = 1 for the speed of light and m = 1/2 for the mass of the particle. 



Let ip he a. state in "H, then it can be uniquely written a.s ip = "V^loop + V^: 



LEAD) 



with 



V'loop ^ -^^([0)-^]) V'lead £ -^^([-^,00)). To simplify the notations, we shall use the 
following convention: ip{x) = ipLoop{x) if x G (0,L), tp{x) = ipLEAoix) if x G {L,oo) and 

1p{0+) = ^LOOP(O), '4'{L-) = ^LOOP(^), = ^LEAd(^)- 

The Hamiltonian H is an unbounded operator and thus cannot be defined on the entire 
Hilbert space "H. We thus have to find the domain D{H) C "H corresponding to the situation 
depicted in Fig. □ such that eU for all ^ G D{H) and {H,D{H)} is self-adjoint. The 
)recise definition of a self-adjoint operator is given in Appendix A and we refer to the paper 



20| for some physical motivations. This section is not crucial to understand the remaining 
of this paper, so the uninterested reader may look at the scattering matrix ([7]) and then go 



directly to the solution (|T9|) . In the paper 



10| . the scattering processes occurring at the 



connection point (between the loop and the lead) are described in terms of the following 
energy-independent (unitary) scattering matrix: 

(-{a + b) 

^/e a b 
y/e b a 



S 



\ 



(7) 



where a = — 2e — l), b = — 2e + 1) and e G (0, |]. Here, e = corresponds to the 
uncoupled situation (which is excluded) and e = | to the maximally coupled one. Note that 
a and b never vanish. To implement such a scattering matrix, it is convenient to work in the 
local reference coordinates {xi G [0, oo) and X2, G [0, L]} associated to the three branches 
exiting from the connection point Xi = X2 = x^ = 0. See Fig. [2j Then, a general scattering 
matrix S{k) at energy k"^ > 0, with k > 0, relates the incoming and outgoing amplitudes at 
the connection point through the following solutions of the stationary Schrodinger equation 
H(f^ = with A; > and £ = 1, 2, 3: 



■ * 

l-f-X2 



i-rX3 



(8) 
(9) 
(10) 
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-|2J], we shall 



where Sij is the Kronecker delta. Inspired from quantum networks studies 
now derive the corresponding boundary conditions at the connection point insuring the 
self-adjointness of the Hamiltonian H. 

X3 



Xi 




FIG. 2. The local reference coordinates. 

Let / and g be two functions in "H such that Hf and Hg are still in H. Then, a simple 
calculation gives 

{f\Hg) = {Hf\g) + BT , (11) 

where (-I-) denotes the usual scalar product in L^([0, oo)) and BT contains the boundary 
terms appearing during the integrations by parts: 

3. r n r n 

(12) 



BT = 5^ [/,(0)^;.(0)-/j (0)^7,(0) 



1j 



/3(0)^73(0)-/2(0)(72(0) 



where /i(xi) = /(L + xi), f2{x2) = f{x2), fsixs) = /(L-xg) and similarly for g^, g2 and 5(3. 
The self-adjointness of H requires that BT vanishes for all f = (/i, /2, /a) and g = {gi,g2, gs) 
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satisfying the same boundary conditions at Xi = X2 = X3 = (see Appendix A for a more 
precise statement). Setting f = g, one thus has the following necessary condition: 



f'(0))c. - (f (0)|f(0))c3 + 2t- [\m\' - l/2(0)P] = 



(13) 



where (■|-)c3 denotes the canonical scalar product in C'^. If || ■ ||c3 = \/(4Tc3, then fll3p 
is satisfied if and only if ||f(0) + 2[f'(0) - if h(0)]||c3 = ||f(0) - ^[f' (0) - if h(0)]||c3, where 
h(0) = (05/2(0), — /3(0)). This is the case if and only if there is a unitary matrix U such 
that 



f(0) 



f (0) - ^|h(0) 



U\ii0)-i 



f'(0)-z|h(0) 



or 



(1 - [/)f (0) + z(l + f/) 



f(0) 



z-h(O) 







(14) 



(15) 



We thus have the following necessary condition for the self-adjointness of {H, D{H)}: there 
must be a unitary matrix U such that the boundary conditions (fT5|) are satisfied for all 
/ G D{H). 

To specify the boundary conditions (|T5|) in terms of the scattering matrix S{k), one 
requires that the scattering states ip^ (which are not in "H) satisfy the relation f|T5l) . Setting 
/ = if^, with i = 1,2 or 3, one easily checks that the relation (|T5|) is verified with U = 
[{k + l)S{k) — {k — 1)1] ■ [{k + 1)1 — {k — l)S{k)]~^, which means that the scattering matrix 
may be written as S{k) = [{k - 1)1 + (k + 1)U] ■ [{k + 1)1 + (k - l)U]~\ Note that the 
knowledge of the scattering matrix S{k) at any energy E = k"^ uniquely determines the 
boundary conditions f|T5|) and that if S does not depend on k, then one simply has U = S. 
Returning to our original reference frame, the scattering matrix (jTj) leads to the following 
boundary conditions: 



V;(0+) - V^(L_) = , 
V^^iO+) - H{L+) = , 
v^^'(L+) + 6[^'(0+)-^'(L_)] = 0. 



(16) 
(17) 
(18) 



Note that these conditions do not depend on the magnetic flux $, that they are contained in 
the model considered in jl6| and that the free boundary conditions ip{0+) = il){LJ) = ^/'(L+) 
and ip'{L_) = ip'{0+) +ip'{L^) correspond to the case e = 4/9. As explained previously, the 
boundary conditions ( 1TB|1 -( IT5|1 are necessary for H to be self-adjoint. In Appendix A, we 
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show that these boundary conditions are also sufficient for H to be self-adjoint. Therefore, 
we deduce that the relations f lTB]) -f lTH]) are in fact the unique boundary conditions associated 
to the scattering matrix ([7j) for which the Hamiltonian ([6]) is self-adjoint. More precisely, 
the Hamiltonian H associated to S is self-adjoint if and only if 

D{H) = {'^ en \ Hip en and - 1^^ are satisfied} (19) 

III. THE PROPAGATOR 

Let us recall that the model we are describing contains three parameters: the loop's length 
L > 0, the magnetic flux $ G M and the coupling parameter e G (0, |]. In what follows, 
we shall suppose that their values are kept constant and that at time t = the initial state 
ipo ^ D{H) of the particle is localized in the loop, i.e. it satisfies ipo{x) = 0, for all x ^ [0, L], 
and is normalized: \tlJo{x)\'^dx = 1. Note that the property tlJo{L+) = ipQ{L^) = together 
with the boundary conditions (fT6|) - (|T8|) imply ipo{0^) = ipo{L^) = and iPq{0^) = iP'q{L^). 
Then, the nonescape probability is given by (t > 0) 

P{i,o.L,^,e,t) = I \ip{x,t)\^dx , (20) 

where the wave function obeys the Schrodinger equation 

d 

i—ip{x,t) = Htp{x,t) , with ip{x,0) = ipo{x) . (21) 

We see that we need to solve the Schrodinger equation (12T!) inside the loop (x G [0,L]) in 
order to compute the nonescape probability fl2U|) . To determine ip{x,t) we shall write the 
initial state 'iPq{ clS cL superposition of generalized eigenstates of H, satisfying (fTBi) - (IT5|l . 
over the spectrum of H, which is crlH) = [0, oo) (see Appendix A). 

Let k > 0, then the general solution of the time-independent Schodinger equation Hipk = 
k'^fk is 

(^fc(x) = [A{k)e"'^ + B{k)e-''^] e^^^X[o,L)(x) 

+ [C(A;)e^'=(^-^) + D(fc)e-^'^(^-^)] X[l,oo)(x) , (22) 

where A{k), B{k), C{k) and D{k) are complex numbers. We shall choose these constants 
such that the boundary conditions (fTB|) - f|TS|l are satisfied and such that the generalized 
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eigenstates {ipki^)} are 5-normalized (as we shall see ip^k{x) = —ipk{x), so only the fk{x) 
with positive k need to be considered): 

POO 

{^k\Vp) = / ^k{x)(Pp{x)dx = 5{k - p) , for all k,p> . (23) 
Jo 



Solving ([IS])-([IH]) and d^S]), one finds |25| 



-i{kL-<S>) _ 1 

A(k) = — ^ , (24) 

i{kL+'S>) _ 1 

1 (26-e)(l + e2**)-2(6-e)e-'('^^-*)-26e'('=^+*) 
^^^^"267! DEN{k) ' ^^^^ 

1 (26-e)(l + e2'*) -2(6-e)e^('=^+*) -26e-^('=^-*) 
^^^^ " "2Vf DEN(^) ' ^^^^ 



where 



DEN(fc) = i^^{a^sm\kL) + [cos(fcL) - cos(<l')]'} 



1/2 



(28) 



In Appendix B, we show that the coefficients A, B, C and D are well defined in M. In 
addition, if cos($) = 1 or cos($) = —1, there are the bound states (p^,(p~ G D{H) of H 
(C = = and A = -B): 

cos($) = l:if+{x) = y^sin (/canx) e**^ X[o,L)(a;) , n = 1, 2, . . . , (29) 

cos($) = -1 : v?„(a;) = W-sin(/i;2„+ix)e*i'' X[o,L)(a;) , n = 0, 1,... (30) 

where /c„ = Note that the eigenvalues En = (x)^' with = 1,2,..., associated to 
the bound states (I29|) - (l30|) are as expected inside the continuous spectrum [0, oo) of H. If 
cos($) = 1, then one has for all m, = 1, 2, . . . and every k > {k E {^^}'^=i included): 

i'Pml'Pt) = ^mn and {<^k\vt) = • (31) 
Similarly, if cos($) = —1, then for all m, n = 0, 1, . . . and k > {k G included): 

(<^ml<^n) = Smn and (<^A,#n) = • (32) 

Since the Hamiltonian H is self-adjoint, it follows from the Spectral Theorem that the 
family of states {v^fc, 5cos($),i ■ , '^cos($),-i ■ Vn.} generates the whole Hilbert space "H. To 
avoid technical complications, we shall assume that the initial state ^pQ can be written 



punctually in terms of these states [this assumption may be checked exphcitly in the cases 
{e = 1/2, $ = 0) and {e = 1/2, $ = n/2)]. 

Assumption (completeness): Let tpQ G D{H) be an initial state localized in the loop. 
Then, we assume that for (almost) every x G [0, oo): 

"^0 n=l 71=0 

Under this assumption and using the orthogonality relations 0231) and fl3ip -( 15^ . one may 
write the solution of the Schrodinger equation (HJ^i as 

-iHt„ 



z^(x,t) = (e-*^Vo)(x)= / K{x,y,t)My)dy , (34) 
where the propagator is 



oo 



K{x,y,t)= e '(pkix)ipk{y)dk 
Jo 

oo 

+5cos{*),i ■^e-*'^'"*<^+(x)<^+(y) 

n=l 

oo 

+^cos{$),-l ■ Yl ix)'fniy) ■ (35) 

n=0 

Looking back at fl2D]) and fl34p one sees that it is sufficient to know the propagator K{x, y, t) 
for x,y E [0, L). In this case, one may write 

K{x, y, t) = Ki{x, y, t) + (5cos(*),i ■ ^2(3;, y, t) + 5cos(<i.) -1 ■ K^ix, y, t) , (36) 

where 

roo 

Ki{x, y, t) = e'*("-^) / \ /,(A;)e-'['='*-'=("-^)l + f^{k)e-'^'''^-''^-+y)] I rfA; , (37) 



2L 

ri=— 00 



2 y ^ 

K3{x,y,t) = ^ ''^^ J2 |e-*['=2Vi<-W(^-y)] _e-*K+i*-W(^'+?/)]| , (39) 



n.=— 00 



with 



T17T 

fcn = , G Z , (40) 

47r {a2 sin2(A;L) + [cos{kL) - cos(<l>)]'} ^ ^ 

^ h +g-2^fcL_2e~ifcLcos($)l , , 

f2{k) = — oT • 42 

87r{a2sin2(A;L) + 62[cos(A;L)-cos(<l>)]'} ^ ^ 
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Clearly, the propagator Ki{x, y, t) describes the decaying aspect of the initial state ipQ, while 
the two other propagators, K2{x,y,t) and K^{x,y,t), describe the stationary aspect of ipQ, 



26| 



i.e. the part of ipo which remains in the loop for all times 
In the case (e = 1/2, $ = 7r/2), one has DEN(A:) = 
propagator exactly by splitting the different terms. We find 



An and one may evaluate the 



K{x,y,t) 



2e 4' 



e 4t 



.{x-y + L)^ 



it -\- le 4* 



(43) 



Also, in the case {e = 1/2, $ = 0), one has fi{k) = l/(47r) and f2{k) = e~*^'^/(47r), so that 



(e=i*=0) 

K{x,y,t) - 



■(x-yy 

e 4* 



(x + y-Ly 

e -it 



+ K2{x,y,t) . 



(44) 



y/lQirit 

In general, we are not able to integrate explicitly Ki{x, y, t), so we shall only discuss its long 
time behaviour (some comments concerning the short-time behaviour of P{ipo, L, $, e, t) can 
be found in Sec. IVII|) . For this, we need to analyse the asymptotic behaviour of an integral 
of the following form: 

/CO 
f{k)e-'^''''-'''^dk , (45) 
-oo 

where / is either fi or /2 and z = {x ± y) G M. Looking at the expressions fHTl) -fjl2 j) . 
(lB5p - (lB6p and (lB8p - (lB9p . it is clear that, for any G M, the function f{k) may be written 
ClS Sb Maclaurin series: 

(46) 



n.=0 



Substituting this expression in ( I45p and exchanging formally the order of integration and 
summation we obtain 



n=0 ' " 



(47) 



We next use the following identity (see Appendix C): 

LfJ 



n-2l 



a j^^il {n-2i)\ (4a)"-^ 



(48) 



where [^J = n/2 if n is even and (n — l)/2 if n is odd, and a,b & C\{0}, with Re(a) > 0. 
The cases 6 = and Re(a) = are obtained as the limiting expressions of the r.h.s. of 
as 6 — )■ and Re(a) — )■ 0"*", respectively. Thus, writing 



e 4t 



E 

m=0 



1 i'^z 



m ^2m 



ml {Aty 



(49) 
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we obtain 



OO OO L 2 J 



22m+n ^! £! (^_2£)! 



1 



n=0 m=0 £=0 

Using the relation (150|) . one can write 

CO OO L'jJ 



^[2(m+n-€)+l]/2 ' 



(50) 



■m.-e-l/2 



1 



22m+n ^! _ 2i)\ t[2(m+n-^)+l]/2 



n=0 m=0 £=0 



(51) 



where the functions /i and /2 are given by fj^Tl) - fj42|) . The first two leading terms are 
K^{x,y,t) = e^^(^-^)| J^[A(0) + ^(0)]^ + {^[/i(0)(a; - + + yf 



'Tl I . 

T 

+ - ^) + + ^)] + ^ [/r(0) + /^'(O)]}^} + O . (52) 

A careful analysis of the expressions fHT]) - fH2|) shows that one has to deal with two cases 
(see also Appendix B): cos($) = 1 and cos($) ^ 1. In particular, one has 

e 



/l(0) = /2(0) 

/i(0) = -/2(0) 



87ra2 



if cos($) = 1 

e 



if cos(<l>) 7^ 1 



(53) 
(54) 



47r62[l - cos($)] 

As we shall see, the fact that /i(0) = /2(0) in the first case, while /i(0) = — /2(0) in the 
second, will lead to a drastic change in the decay properties of the nonescape probability. 
If cos(<l>) = 1, then 



K{x,y,t) = K2{x,y,t) 



(55) 



+ 4 



+ - L{x + y) + Y {- 



1 



o 



^5/2 



If cos(<l>) 7^ 1, then 



K{x,y,t) = 5cos(<i.) -1 K-i{x,y,t) 



le e L 



"Kl 



IP 



(56) 



{L - [1 - cos($) + i sin($)]x} {L - [1 - cos($) - i sin($)]?/} 1 



O 



^5/2 



[1 - C0S($)]2 

As one may easily check, the expression (135]) with {e = 1/2,$ = 0) and the expression fl56p 
with (e = 1/2,$ = n/2) coincide with the leading terms of (jS]) and (HSjl . respectively. 
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IV. GENERAL RESULTS 



From the relations (I55p - fl56p . it follows that 

r T.n .^ / Poo(V^o, I^, if I cos(<l>)| = 1 

otherwise 



where 



Poo(^/'o,^,$) = / dx 



L 







L 



dy [4os(<i.),i K2{x, y, 0) + 5cos(*) -i ^3(3;, 0)] ?/'o(2/) 



2 



(58) 



If cos($) = 1 [or cos($) = —1], then the initial state ipQ may be written as a sum of two terms 
[see 0331) ]. one associated to the scattering states (fk and the other to the bound states (p^ [or 
ip~]. The part associated to the bound states will lead to a non-zero value of -Poo(V^O) L, $). 
Note that Poo{ipo, L,^) does not depend on e. From the expressions fl55|) - fl56|) . one sees 
that the large-time behaviour of P{ipOi L, $, e, t) depends drastically on the value of cos($). 
Nevertheless, setting Poo(V^05 L,^) = whenever [ cos($)[ 7^ 1, one may write 

P(^o, L, $, e, t) = PM,, $) + ^^^^"^f'"^'^^ . (59) 

k=l 

To obtain explicit expressions for the coefficients Cki'ipo-, L,(^,e), it is convenient to write 
the expression fl3T]) in the following form: 

^i(-'^'^) = E§fe|- (60) 

£=0 

Then, the decaying part of the wave function reads 

Ki{x, y, t)ipo{y)dy = ^ .(^2e+i)/2 ' (6^) 



where 

"L 



Qe{x)= I Ge{x,y)My)dy . (62) 
Jo 

Hence, 

P{i:o,L,<l>,e,t)-PUi'o,L,<i>)= f \i^^...,{x,t)\Hx = ^^^(^^0' ^) , (53) 

"^0 fc=i 

where 

00 

Cfc(^0,i^,$,£) = 5^ {Qm\Qi) . (64) 

£, m = 
£ + m + 1 = fc 
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Note that if Cii^ipo, L, $, e) = 0, then Qq = and thus 020^0: L, $, e) = 0. More generally, if 
Ci{ilJo,L,<!>,e) = C2{tpo,L,<!>,e) = ■■■ = C2i+i{^po, L,<!>,e) = 0, then C2i+2{4^o, L,<l>,e) = 0. 
Thus, at large times the nonescape probability P{ipo, L,^,e,t) ~ Cs{ipo, L,^,e)/t^ , with 
6 an odd number. Let us use the expressions ( I55l) - (l56l) to write the leading coefficients 
Ck{ipo, L,^,e). It turns out that they can be nicely written in terms of the derivatives of 
the Fourier transform of the initial state: 

= r e-^'yy-My)dy • (65) 

V 2lT Jo 

There are two cases: 



;i) If cos($) = 1, then 



Ci{ipo,L,<^,e) 



^0 



If Ci{ipo, L, $, e) = 0, then C2{i'o, L, $, e) = and the next leading term is given by 



128 a4 



(2) If cos(<l>) ^ 1, then Ci{ipo, L,^,e) = C2{ipo, L,^,e) = and 



C3{ilJo,L,^,e) 



96 b^[l - cos(<l))]4 



(66) 



(67) 



(68) 



{sin(<l>) + i[l - cos(<l>)]}^o 



V'o-dependence: We see that only the momentum k = ^/L plays a role in the decaying 
properties of Pltpo, L,^,e,t). When $ = 0, this coincides with the well-known fact that 
only the zero momentum is involved 3|. Note that Ci{ipo, L, $, e) 7^ only if cos($) = 1 and 
ijjo{^/L) 7^ 0. As an illustration, let us discuss the case $ = 0. In this case, the bound states 
(BUj) form a basis in L'^{[0, L]) of all odd functions with respect to L/2 [i.e. ip{x) = —ip{L — x) 
for all X G [0,L]]. Thus, if $ = and the initial state ipo ^ D{H) is an odd function with 
respect to L/2, then P(^/;o, L, $, e, t) = Poo{'ipo, L,(^) = 1 for all times. On the contrary, if 
$ = and ipo G D{H) is an even function with respect to L/2 [i.e. iPq{x) = ipo{L — x) for 
all X e [0,L]], then Poo(^/'o, ^, = and P{4jo, L,<l>,e,t) = Ci{tpo, L,<l>,e)/t + 0{l/t'^) at 
large times, with Ci{i/jo,L,^,e) 7^ since Jq ipo{y)dy 7^ 0. More generally, if $ = and 
i^o = ^S^'^ + i^r"^ then Pi^Po, L, <D, e, t) = P^i^Po, L, $) + Ci(^Ao, L, <D, e)/t + 0{l/t^), where 
Poo(^o,i>,$) = Iodx\ jI'dyK2{x,y,0W{y)\' and Ci(^/;o, L, <D, e) = e^L/{Sa'') |V5§^^"(0)|2. 
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2y-dependence: Let iIjq{x,X) and K{x,y,t,X) be the initial state and the propagator 
corresponding to a loop of length XL, with A > 0. Then, from the relations fl72]) and fl5^ . 
one deduces that 



From the relation (159|) . one easily deduces that Poo is independent of L, so from now on we 
shall write Pooi'ipo, and that the coefficient Ck is an homogeneous function of degree 2k, 
i.e. Ck{i>o, XL, $, e) = X^'' Ck{ipo, L, $, e) for all A > 0, and thus Ck{ipo, L, $, e) oc L^^. 

^-dependence: Looking at the relations fl66l) - fl68l) . one sees that in general it is difficult 
to predict the ^-dependence of P{ipQ, L,^,e,t) and we shall discuss this matter in detail 
in the next section. Nevertheless, some general symmetries can be found. Indeed, the 
propagators Kj(x, y, t), with j = 1, 2, 3, become e^^'^^'^^ Kj{x, y, t) under the transformation 
$ $ + 271. This implies that P^ipo, L,^,e,t) is invariant under this transformation on 

■ 27V 

the proviso that one replaces the initial state ipo{y) by e~''~'^ipo{y). Furthermore, a simple 
calculation shows that P{ipo, L, $, e, t) is invariant under the transformation $ i— t- — $ if the 
initial state satisfies ipoi^) = ~V^o(-^ ~ [as in f l72]) ] or ipoi^) = i'oi.L — x), for all x G [0, L\. 

^-dependence: At large times, one sees from fl^B]) - fl^5]) that P{ipQ, L, $, e, t) — Poo^ipQ, $) 
behaves as e'^/a^ if cos($) = 1 and e'^/b^ otherwise. Looking at Fig. |3l one deduces that 
at large times P{ipo, L,^,e,t) — Pooi'ipo,^) decreases with e if cos($) = 1 and increases 
otherwise. In other words, in the long run the particle escapes faster as one increases e if 
cos($) = 1 and escapes slower otherwise. This slowing effect when cos($) 7^ 1 at large times 
is rather peculiar since in the short time regime, some numerical results show (see Fig. [TUj) 
that the particle may escape faster as e increases. 



\/Xipo{Xx,X) =^o{x,l) , 
X K{Xx, Xy, XH, A) = K{x, y, t, 1) . 



(69) 
(70) 



From these scaling laws, one finds that 




(71) 
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0.4 0.5 



FIG. 3. Log-lin plot of e^ja^ (dashed blue) and l}f (red) for e e [0,0.5] 



V. THE BOUND STATES 



In order to apply our general results, let us consider the following situation: the system 
is prepared such that at time t = no magnetic flux crosses the loop and the particle is in 
a bound state V9„ G D{H), where 



|sm 



27m 



L 



X X[o,L){x) , n = 1,2. . . . 



(72) 



Since the effect of the length's loop and of the coupling parameter was thoroughly discussed 
in Sec. IIV| let us fix L and e and discuss the ^-dependence of P{(pn, L, e, t). If at all later 
times there is still no magnetic flux inside the loop ($ = 0), then all the bound states (pn 
will lead naturally to P{ipn, L,^,e,t) = 1 at all times, i.e. the particle will remain forever 
in the loop. On the contrary, if some constant magnetic field is applied ($ ^ 0), then the 
states ipn will no longer be bound states of the Hamiltonian and consequently the particle 
may escape from the loop. Using 0581) . one finds the following asymptotic values: 



lim P(y?„, L, e, t) = P^{^n. $) 



1 if $ = 
if |cos($)| 7^ 1 
1/4 if |$| = 2mx 
^1/2 otherwise 



(73) 



As explained in Sec. HV] P{(pn, L, e, t) is invariant under $ i— )■ — $ since ipn are odd 
functions with respect to L/2, so it is sufficient to discuss the case $ > 0. If $ = 0, then 
P{(Pn, L,^,e,t) = 1 at all times. If $ G (0,27r), then by using ([73]) and ([SS]) one finds 
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P(^„,, L, $, £, t) = 5$, J2 + C3(</^n, L, $, + C(l/t4), where 

2nV$2L67r[2 + cos(<l>)] 



C3(y?„,L,$,e) 



364[1 - cos($)]2($ - 2n7r)4($ + 2n7r)4 ' 



(74) 



A simple analysis of (17^ reveals that C3 ((/)„, L, $, e) is always positive and diverges at the 
excluded values cos($) = 1. Furthermore, we see in Fig. IHthat C^{ipn^ L, $, e) first decreases 
and then increases with $. Note also that C^i^Lpn, L, $, e) decreases with n, showing that in 
the long run more energetic bound states escape faster from the loop. If $ = 27r or more 
generally if cos(<l') = 1, then by using f l66|) - fl67|) one finds P{ipn, L,^,e,t) = Poo{ipn,^) + 
Ciiipn, L, $, e)/t + C2{ipn, L, $, e)/t' + Cs{ipn, L, $, e)/t^ + 0(l/t^), where 



Ci((^„,L,$,e) 



L2e7(327ra^) if |$| = 2mT 
otherwise 



(75) 



If Ci ((/:>„, L, $, e) = (i.e. |$| 7^ 2n7r), then C2(v?n, L, $, e) = and 

C3((^„,L,$,e) 



(76) 



2a4($ - 2n7r)4($ + 2r27r)4 ' 

Note that, as in (I71|) . the coefficient C3(v9„, L, $, 5) in (1761) decreases with n. Note also that 
C3{ipn, L, $, 5) decreases with $, as one may easily see by looking at the red dots in Fig. |5l 
Finally, we see in Fig. |5]that by increasing further the magnetic flux $ we obtain a series of 
similar patterns in each range [2£'it, 2{£ + l)vr], £ = 0, 1, 2, . . . 



e 

I 

U 




0.001 



FIG. 4. Log-lin plot of Csi^pn, L, e) in ([71]) for n = 1, L = 1, e = 1/2 and $ G (0, 27r). 
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20 40 60 80 100 



FIG. 5. Log-lin plot of C3(v3„, L, e) for n = 1, L = 1, e = 1/2 and <P G [0,100]. When 
$ = 0, C3{ipn=i,L = l,^>,e = 1/2) = 0. When $ = 2tt, Ci((^„=i,L = l,$,e = 1/2) / so 
C^{ipn=i-,L = l,<I>,e = 1/2) is not plotted. When $ = 47r,67r, . . . , we use ([75]) to obtain the red 
dots. All the other values are obtained with ()74p . 

To obtain a global picture in time, we have computed P{(pn, L, $, e, t) numerically. For 
this, we wrote f l57|l as 



Ki{x,y,t) 



G{k;x,y,t)dk = 2 / C'ik; x,y,t)dk 
Jo 



where 



G'{k; X, y, t) = -[G{k; x, y, t) + G{-k; x, y, t)] 



(77) 



(78) 



so that 



OO POO 



^0 



L pL pL 



^0 ^0 



G^{k; X, y, t)G^{p; x, z, t)t(jo{y)'ilJo{z)dxdydz 



dkdp . 
(79) 

Then, we integrated analytically the expression inside the square brackets and then com- 
puted numerically (with Matlab) the so obtained expression with respect to k and p over 
the finite domain [0, c] x [0,c]. The constant c was chosen such that at very short times 
P{(Pn, L, $, e, t) was close to one and no appreciable changes occurred in P{(pn, L, $, e, t) by 
increasing c further. We also checked that the curve for {e = 1/2, $ = 7r/2), for which we 
know the propagator exactly [see f H3|) ]. was accurately reproduced. In all the curves having 
reached the asymptotic regime, the long time values are in very good agreement with the 
ones given by the relations ( 171|1 -( 17^ . 
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FIG. 6. Log-log plot of P{ipn,L,^,e,t) for L = 1, e = 1/2, $ = tt/2 and n = 1 (black circle), 
n = 2 (red square) and n = 3 (blue triangle). 

^-dependence: The ^-dependence of P{y:>n, L, $, e, t) is represented in Figs. In the 
curves having reached the asymptotic regime, one sees that P{t) ~ C{ipn, L, $, e) /t^ at large 
times, where the constant C(</9„, L, $, e) decreases with $ when $ G (0,2.33) (see Fig. [7]) 
and increases when $ G (2.33, 27r) (see Fig. |8]), as expected from Fig. HI Note however that 
there may be several crossing among the curves. Note that similar features were obtained 



in a circular dielectric cavity containing classical waves [19|. In Fig. |9l one sees that $ = 27r 
leads to P{(pn, L,^,e,t) ~ C{(pn, L,^,e)/t as expected from fl75]) . while $ = 7r,37r,47r lead 
to P{(pn, L,^,e,t) ~ C{ipn, L,^,e)/t^ as expected from ([71]) and (17^ . 

^-dependence: In Fig. [TOl one sees that P{(pn, L,^,e,t) ~ C{(pn, L,^,e)/t^ at large 
times and that the constant C((y9„,, L, $, e) increases with e. On the other hand, note that 
P{(fni L, $, e, t) decreases with e at short times and that the curves cross each other. 

Exponential decay: Plotting the curves in Figs. [6]-[lQ] on a log-lin plot (not shown) re- 
veals that on some intermediate time scale (not very small and not very large), the nonescape 



probability decays approximately as an exponential, in agreement with 



i 
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10"^ 10"^ 10° 10^ 10^ 

t 



FIG. 7. Log-log plot of P{ipn,L,^,e,t) for ?i = 1, L = 1, e = 1/2 and = 0.1 (green triangle), 
^ = 0.3 (yellow diamond), <I> = 0.5 (purple circle), $ = 0.7 (pink lozenge), <^ = 1 (red square) and 
$ = 2 (blue cross). 




-2-10 1 2 

10 10 10 10 10 

t 



FIG. 8. Log-log plot of P{ipn,L,^,e,t) for n = 1, L = 1 e = 1/2 and <I> = 3 (red circle), <I> = 4 
(green square), $ = 5 (black triangle) and $ = 6 (blue cross). 
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FIG. 9. Log-log plot of P{(pn,L,^,e,t) - Poo{fn,^) for n = 1, L = 1, e = 1/2 and $ = vr (black 
circle), $ = 27r (blue square), <I> = Svr (red triangle) and <1> = 47r (green diamond). 




VI. THE INTERACTING CASE 



So far we have solely discussed the free evolution. To have an idea about the interacting 
evolution, let us consider the case $ = and e = 4/9, where the Hamiltonian of the particle 
is given by 



H = Ho + V{Q) 
21 



(80) 



Here Hq = — ^ is the free Hamiltonian and V is the potential which for simphcity is 
supposed to be bounded, so that the Hamiltonian H is self-adjoint in the domain D[H) 
specified by the conditions (ITB]) - (IT5]) with e = 4/9. 

Obviously, in this case the time- dependent Schrodinger equation (l2Ti) cannot be solved 
analytically in general and it also seems extremely difficult to compute analytically the large 
time behaviour of the nonescape probability, so we shall consider a fully numerical approach. 
For this, let us consider a long but finite lead of length i > L and set ip{L + £) = 0. The 
assumption i > L is natural in our case and also simplifies the remaining discussion. One 
can then check that H is self-adjoint in 

D{H) = {ij e n \ Hij e n and - are satisfied} (81) 

where 

^(0+) = ij{L_) = i,{L+) , (82) 

V^'(L_)=V^'(0+) + ^A'(L+) , (83) 

i/j{L + £) = 0. (84) 

Clearly, for our numerical results to remain accurate for a long time the unitarity of the 
time-evolution is crucial. A convenient way to implement numerically the boundary condi- 
tions fl82l) - fl84l) . while preserving unitarity, is to use the pseudo-spectral method presented 
in Appendix D. This method requires the knowledge of the spectrum and bound states of 
the free Hamiltonian Hq, which are worked out in Appendix E. We believe that this section 
and the Appendices D and E present, in particular, an interesting and highly non-trivial 
application of the pseudo-spectral method. 

Let us consider the following initial state: 

2 /27r 




M^) = V"^ / X[o,L]i^) ■ (85) 

It is easy to check that tpQ is normalized (HV^oll = 1) s-nd belongs to D{H). Note that ipQ 
is an even function in L^([0,L]) and thus it is orthogonal to all bound states of the loop. 
Therefore, in the infinite case i = oo with no potential (V^ = 0), the nonescape probability 
will vanish at large times and the asymptotic decay can be computed by using fl66|) : 



Pfree(^0, L,i = 00,t) = — - + Oi-\ . (86) 
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Numerically, the infinite situation {i = oo) is obtained by using the expression (17^ and 
gives the dark blue cross curve in Fig. [TT] Its asymptotic values are in very good agreement 
with formula fl86|) . On the other hand, the finite case {£ < oo) is obtained by using the 
pseudo-spectral method and leads to the light blue curve in Fig. [TTl As one clearly sees, the 
two methods are in very good agreement. 





10°' 








10-' ■ 








10-^- 


4-' 




Q_ 


10-'- 








10" ■ 








10-=- 

1 




10 




-3 



10 



-2 



10 



-1 



10 



10 



FIG. 11. The nonescape probability P{^o, L, i, X,t) for a soft-core coulomb potential by using the 
pseudo-spectral method with L = 1, £ = 500, Ax = 1/80, At = 10"'^. From top to bottom: 
A = —0.1 (green), A = (light blue), A = 0.1 (orange) and A = 1 (red). The dark blue cross curve 
(close to the A = curve) corresponds to the free case {V = 0) with i = oo and is obtained by 
using the expression ([79]) . 



Imagine now that a fictive charged particle is fixed at the point x = L/2, i.e. at the 
leftmost part of the loop, and interacts with our genuine charged particle via a soft-core 
coulomb force. Then, a simple calculation leads to the following potential {0 < x < L + £): 



Vx(x) 



A 



d{x) 



— cos (tt — ^) if X < L 
2R + X- L iix> L 



R 



When A > 0, the force is repulsive, while A < corresponds to an attractive force. The 
numerical results are presented in Fig. [TT] As expected, we see that Pi'ipo, L, i, A, t) decays 
faster than the free case if A > and slower if A < 0. Note also that globally the curves look 
similar to the ones obtained in the free cases presented in the previous section. In particular, 
plotting these curves on a log-lin plot (not shown) reveals the approximate exponential decay 
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on some intermediate time scale. Although our numerical results do not show any asymptotic 
power law decay when A 7^ 0, such a regime may well be present at larger times. To reach 
larger times, however, one has to increase the length of the lead I (because the particle gets 
reflected at the end of the lead and comes back into the loop causing a bump in the curves 
after some time) but this increases the computation time very quickly. 

VII. CONCLUDING REMARKS 

Although more energetic bound states typically escape faster from the loop, we found 
that different bound states may actually decay with different power laws and consequently 
it may be possible to fllter one bound state from a coherent superposition by applying the 
appropriate magnetic flux $ through the loop. Indeed, suppose that the system is prepared 
such that at time t = no magnetic flux crosses the loop and the particle is in the state 



satisfying \\iPq\\ = 1, where ipn are the bound states (172|) . Then, if one wants to filter the 
bound state ipm, one may apply a magnetic flux $ = 2rmT through the loop for a sufficiently 
long time, so that the state (fm decays as 1/t while all the others decay asl/t^ [see (!75|) -( I76|) ]. 

The short-time behaviour of P^ipQ, L,^,e,t) is also interesting and can be conveniently 
written as P{ipQ, L,^,e,t) = 1 — | |x[L,oo)(<5)e~*''^Vo| By computing the propagator 
K{x, y, t) for X > L and y < L, one then founds that P{ipo, L, $, e, t) = l—D2{ipo, L, $, + 
0{t^) as t — 7- 0. In general, it seems that the constant D2{ipo^ L,^,s) cannot be expressed 
in a simple form and thus its behaviour would require a refined numerical investigation. 

To make more apparent the exponential decay of the nonescape probability at inter- 
mediate times, while keeping the power law decay explicit at large times, it would be 
interesting to write P{'ipo, L,^,e,t) as Poo{ipo,^) + A{ipo, L,^,e) exp[—a{'ilJo, L,^,e)t] + 
Y^'^=iBk{ipo,L,^,e)/t'' and to determine the unknown variables A, a and Bk- One may 
also analyse the transition times between the exponential regime and the power law regime 
with respect to the parameters e and $, and also A if there is such a transition. 



00 




(88) 



n=l 
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APPENDIX A: PROPERTIES OF THE HAMILTONIAN H 



Here we show that the Hamiltonian H is self-adjoint in D{H) given in f|T9i) . First of all, 
it is clear that D{H) is dense in "H, so that the adjoint H* of H is well defined. Let us 
recall that an operator {H,D{H)] is self-adjoint if D{H*) = D{H) and H*il) = Hip for all 
ip G D{H), where a vector / G "H belongs to D{H*) if there is a vector /* G "H such that 
{f\Hg) = {f*\g) for all g G D{H), and for each / G D{H*) one sets H*f = f*. 

Let / G "H be a vector such that Hf G "H and let g G D{H). Then, 

{f\Hg) = {Hf\g) + BT , (Al) 



where 



BT 



2i-g{Q^ 



) /(^-)-/(o+) 



g{Q^ 



{v^/'(L+) + 6[/'(0h.)-/'(L_)]} 



v/F/(0+) - hf{L, 



Here we have used the fact that g G D{H). We see that BT vanishes for all g G D{H) if and 
only if / also satisfies the boundary conditions (fTB]) - ffTS]) . This shows that D{H*) = D{H) 
and H*f = Hf for all / G D{H), i.e. H is self-adjoint. 

Let us now show that the spectrum (t{H) of H is [0, oo). Since H is self-adjoint, it 
follows that o-{H) C M. Writing H = P|, with P$ = — — f X[o,l](Q); one easily shows 
that {(p\Hif) = ||P$(y9|p > for all if G D{H). This means that if is a positive operator 
and thus cr(if ) c [0, oo). Let A; > be fixed, then a (particular) solution of Hip = k'^ip is 

y,(x) = e^'=(--^)x[L,oo)(x) . (A2) 



Take ip G Co~(M) satisfying ip{^) = ^^'(0) = and \ij{y)\^dy = 1, and set 



n 



(p{x) , where a: > and n 



1,2, 



(A3) 
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It is easy to check that ipn G D{H) and = 1 for all n, and that {{Hipn — k'^ipnW — as 
n — 7- oo. This shows that E = k"^ > belongs to the spectrum of H (see Prop. 4.20 in 

28|). Since > is arbitrary and the spectrum of if is a closed set, one can 



or Th. 2.2.1 in 



conclude that cr{H) = [0, oo). 



APPENDIX B: THE COEFFICIENTS A, B, C AND D 

The denominator DEN(A;) of the coefficients fl2^ - fl27|) vanish if and only if 

sin2(fcL) + 62[cos(A;L) - cos{^)f = . (Bl) 

Since a and b are always non-zero (remember that e = is excluded) the above equality is 
satisfied if and only if 

sin(A;L) = and cos(fcL) = cos($) . (B2) 

The first relation is verified if only if kL = nn, for some n G Z, and consequently the second 
relation reads cos(<I>) = (—1)". This shows that the coefficients A, B, C and D may contain 
singularities only if cos(<l>) = 1 or cos($) = —1. We shall now discuss these two cases 
separately. 



(i) Case cos($) = 1 or $ = 2mi^ with n G Z: Using trigonometric relations, one may 
write 



e a/1 — cos{kL) + i sign[sin(fcL)] ^1 + cos{kL) 



A(k) = -B(k) = -J—^ ^ ^ ^ _v 

VStt {a2[l + cos(A;L)] +62[i _cos(A;L)]}'/' 

Since the denominator of this expression is always non-zero, it follows that A and B are well 

defined in M. Next, one may write 

1 



C{k) = -D{k) = -^[2(& - e)A{k) - 2bB{k)] , (B4) 
showing that C and D are also well defined in M. Furthermore, one has 

^'^^^ " 47r{a2[l + cos(A;L)] + ^^[i _ cos(fcL)]} ' ^^^^ 

^'^^^ " 47r{a2[l + cos(A;L)] + b^[l - cos{kL)]} ' ^^^^ 
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(ii) Case cos($) = — 1 or $ = (2n + l)-7r, with n ^"L: In this case, the relation (IB4p still 
holds with 



e a/1 + cos{kL) — i sign[sin(A;L)]A/l — cos{kL) 



A(k) = -B(k) = -J — ^ ^ — ^ — ^ — '-^ — . (B7) 

V Svr {a2[l -cos(A;L)] +62[l + cos(A;L)]}'/' 

Therefore, the coefficients are also well defined in this case and we have 

•^'^^^ " 47r{a2[l - cos(A;L)] + b^[l + cos{kL)]} ' 

^'^^^ " 47r{a2[l - cos(A;L)] + b^[l + cos(A;L)]} ' ^^^^ 
APPENDIX C: THE IDENTITY (f48]> 
Here we show the identity (jUD- From the equality 



OO /_ 2 



e-'='+^^-rf/t = J-ei^ (CI) 
V a 



one can write 



One then concludes by using the following formula for the Hermite polynomial: 

I - 1 

£=0 ^ ^ 

APPENDIX D: THE PSEUDO-SPECTRAL METHOD 

In Appendix E, we show that the free Hamiltonian Ho has a pure point spectrum {En}^=i 
satisfying < Ei < E2 < ■ ■ ■ (so each eigenvalue is positive and has multiplicity 1), and 
that one can solve (explicitely) the free stationary Schrodindger equation 

with ipn e D{Ho) = D{H) and = [J^^^ \^n{x)\^dxY^^ = 1. Let Vo e D{H), x G 

[0, L + i] and At > 0, then 

V'(x,At) = (e-'^^Vo)(a:) • (D2) 
Let us consider the symmetric decomposition known as Strang splitting {H = Hq + V): 



-iHAt 
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Since {y:>n}n>i forms an orthonormal basis of the Hilbert space "H, one can write 

oo 

e-^^*/Vo = $^(¥.„|e-^^^*/Vo)¥^n . (D4) 

n=l 

Substituting dDS]) and in gives 

ilj{x, At) = J2 e-^^-^*e-^^(^)^*/Vn(a:) / e-'''^y^^'/'^r.{y)My)dy + 0{At') . (D5) 

n=l 

Let $ G "H, and let us define a generalized Fourier transform J-" as (n = 1,2,...) 

l>(n) = 7{^){n) = / ^n{ymy)dy . (D6) 
Since the family {v9„}„>i forms an orthonormal basis of "H, we have (x G [0, L + £]): 

oo 

^{x) = 7-\^){x) = <^m(a;)<l(m) . (D7) 

m=l 

One may then rewrite f IDSp as 

iIj{x, t + At) = e-^^W^*/2jr-i [g-^s.At . jr ^g-iy(.)At/2^^.^ ^^^j j ^^^j ^ 0(At=^) . (D8) 

This relation is a common way of presenting the pseudo-spectral method. Numerically, the 
Fourier transform and its inverse are then computed by using a very efficient algorithm 
known as the Fast Fourier Transform (FFT). As far as we know, there is however no FFT 
using the eigenvectors {(pn}n>i of our model, so we shall compute ip{x,t + At) directly. 

Let Ax = {L + i)/N, with G N, be a small length difference, Xk = k Ax the position 
of the A;-th site, with k = 0, N, tm = rn At, with m = 0, 1, . . . , and 0) = ipQ. Then, 
we will use the following relation in the numerical computation: 

i;{x,, Wi) = E e-^^"^*e-^^(^^)^*/Vn(x.) / e-^^(^)^*/Vn(l/)^(y, t„.)dy + 0{At') . 

n=l -^0 

(D9) 

In the simulations, we consider only the first 2500 terms in the infinite sum occurring in 
( ]D9p . compute the integral by using the composite Simpson rule and check that probability 
is conserved at all times. 

APPENDIX E: THE FREE HAMILTONIAN Hq 

The free Hamiltonian Hq is self-adjoint and thus its spectrum cr{Ho) C M. Writing 
Ho = P^, with P = one finds that {ip\Hotp) = ||P<^|p > for all ip G D{H) and thus 
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(7{Hq) C [0, oo). The general solution of -ffoV^n = Enipn, with En = and /c„ > 0, is 

^n{x) = [Ane'""'' + B^e-''-^] X[o,l)(x) + [C„e^^"^ + D„e"^^"^] xlL,L+i]{x) , (El) 

where X/ is the characteristic function of the interval /: Xi{x) = ^ x & I and Xiix) = 
otherwise. Solving fl82l) - fl8^ . one gets two types of solutions: First, we have the bound 
states of the loop (C„ = Dn = and An = —Bn): 

/2 2TTm 
-sin(Cx)x[o,L)(x) , C^ = -^, m = l,2,... (E2) 

Second, assuming that kn 7^ all m = 1, 2, . . . , we have 

Bn = -Y^^^n , (E3) 
Q _ Ap-iknL _|_ p-2iknL 

= 2(1 -eXx) ■ 

Dn = - Jf . 1 , (E5) 



2(1 -e- 




3 _ 4(.iknL 


_j_ g2ifci 


2(1 -e 


-iknL^ 


(2s + l)7r 
2i 


s = C 



Explicitly, we thus have 



= I ^%r.L {2 [sin(A;„x) - sin(A:„(x - L))] X[o,l){^) (E8) 
+ [3 sm{knx) - 4 sin{kn{x - L)) + sin(/c„,(x - 2L))] xlL,L+i] (x) } , 

where the coefficient An is determined by normalization: J^~^^ \(pn{x)\'^dx = 1. After some 
algebra, we found 

sin (knL) [1-2 sin^ (kni)] 



An = <2L+ £[5-3 cos (knL)] + 2 



^ sin (fc„£) cos (kni) [5 cos (fc^L) - 3] | ^^^^ 



It remains to solve (jETj). As far as we can see, this relation cannot be solved explicitly. 
Since we shall sum up over the energies En = /c^, we actually only need to find out the non- 
negative solutions of (IE7p . Let Iq = [0, ^] and I = [—fi,-^], and let us look for solutions of 
the following form: 

kn = nj + Sn , (ElO) 
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where n = 0, 1, 2, 
for all m = 0, 1, 2, 



, Sq G /q, Sn ^ I and such that /c„ 7^ 211171/ L and /c„ 7^ (2m + l)7r/(2£) 
. . Then, the relation flE7D reads 



2 tan(n7r + 



sin(n^ + 6nL) 
1 - cos(n^ + 6nL) 



(Ell) 



Since tan(n7r + = tan((5„£), one finds 



- arctan 



sm m 



ttL 



+ SnL) 



2[1 - cos(n^ + 



(E12) 



where we consider the principal value of the arc-tangent function: arctan('i9) G (— f,f)- Note 
that cos(n^ + ^n-L) 7^ 1, since /c„ 7^ 2'Km/L for all m = 0, 1, 2, ... , so the above expression 
is well defined. Nevertheless, it will be convenient to attribute values to this expression on 
these exceptional points. 

The condition cos{n^ + SnL) = 1 reads + 6nL = m27r, with m = 0, 1, 2, ... , or 
Sn = m2TX / L — n-nji. Note that 5o £ -^0 if and only if m = 0, since £ > L, and thus 5o = 0. 
For 5„ to belong to /, with n > 1, it is necessary that |4m£/L — 2n\ < 1. Let us assume 
that i/L is an integer (as in the simulation), then only the condition \4mi/L — 2n\ = or 
5„ = might be satisfied, and thus any exceptional point 5„ must be at the centre of /. Let 
us introduce the set S = {2ml/ L \ m = 1, 2, . . . }, which is the collection of n for which the 
condition \4mi/L — 2n\ = 0, with m = 1, 2, . . . , is satisfied. If n G S, then we shall cut the 
interval I into two parts: / = U where = [-^,0] and = [0, ^]. 

Let us define Fq : Jq — ;> Iq as 

sin((5L) 



Fo{S) = 7 arctan 

i l2[l — cos(oL)J 

and set -Fo(O) = tt /{2t). If ri > 1 and ri ^ S, then we define F„ :/—)■/ by 



(E13) 



Fn{5) 



1 



arctan 



sm ra- 



5L) 



2[1 - cos(ri^ + 5L)] 



(E14) 



On the other hand, if n G S, we shall consider the functions : ^ and : — 
given by 



1 



F (^) = £ arctan 
F^ {6) = - arctan 



sin(5L) 



_2[l-cos((5L)]_ 

sin(5L) 
_2[l-cos(5L)]_ 



(E15) 
(E16) 
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and set F^(0) = -vr/(2£) and F^(0) = 7r/(2£). 

We shall show that Fq, Fn, F^ and F^ are contracting maps and thus, by the Fixed Point 
Theorem, there exist unique points 6q, 5*, 6^* and 6^* such that Fo{6q) = 6q, Fn{S^) = 5*, 
F^{6^*) = 6^* and = 6^*. Furthermore, these fixed points may be reached by 

iterations, e.g. F^(0), the m-th iterate of F„ at point 0, converges to 5* as m — )■ oo. This 
shows in particular that we have found all the solutions of flE12p and consequently of flE7l) . 
that Ho has pure point spectrum and that each eigenvalue has multiplicity 1. Let us now 
show that Fn contracts, the proofs for Fq, F^ and F^ are similar. So let X,Y E I with 
X ^Y. Then, 



\F^{X)-Fn{Y)\<C\X-Y\ 



(E17) 



where 



^ IFJX) - FJY)\ \Fn(X) - FJY)\ 

C= ' ^' < sup ' ^' <sup 



dF 



dX 



<^<i. 



(E18) 



The last two inequalities follow from the assumption L < £ and the explicit computation of 

F'ixy. 



dF„ 



dX 



{X) 



2L 



i[5-3cos{n^ + LX)] 



(E19) 



It remains to check that the fixed point 5* of -F„ gives rise to an allowed value of 
i.e. satisfying kn ^ {2'n-m/L, (2m + l)7r/(2^)}J^^Q. These conditions state that 5* should not 
be a boundary value of J, i.e. 7^ 7r/(2£). This can be easily checked. 

In the simulations, we iterate the map Fq, F^, F^ and F^ until the condition f l84p with 



ypsiobai gjygj^ jj^ (|E8|) is accurately satisfied: we required that \lp^°^^\L + i)\ < 10 



Since the contracting constant C is very small (< L/i), this condition is satisfied after a only 
few iterations. In summary, we have 



kQ = Sq , with Sq ^ Fqo ■ ■ ■ o Fo(e) 



77/ ^ Xj ! /c^ 

neT^-.k]:^ 

neT.:k^ 



vr 



n- 



vr 



5: , with 5l^FnO---o F„(0) 



n- + 5^* , with 5 

TT 



5^* , with 5^* ^F^o...o F^(e) 



(E20) 
(E21) 
(E22) 
(E23) 
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where e is a small positive number [e < 7r/(2£)]. The states can then be written as 

if n ^ S : = , 



with k„ 



(E24) 



if n G S : < 



loop 



with kn = f 



^ ^global^ 



with k„ 
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